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Abstract 

We study the phase structures of M = 4 U(N) super Yang-Mills theories on 
M x S'^fZk with large N. The theory has many vacua labelled by the holonomy 
matrix along the non-trivial cycle on 5 3 /Zfc, and for the fermions the periodic and 
the anti-periodic boundary conditions can be assigned along the cycle. We compute 
the partition functions of the orbifold theories and observe that phase transitions 
occur even in the zero 't Hooft coupling limit. With the periodic boundary condition, 
the vacua of the gauge theory are dual to various arrangements of k NS5-branes. 
With the anti-periodic boundary condition, transitions between the vacua are dual 
to localized tachyon condensations. In particular, the mass of a deformed geometry 
is compared with the Casimir energy for the dual vacuum. We also obtain an index 
for the supersymmetric orbifold theory. 
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1 Introduction 

Recently the thermodynamics of large N gauge theories on compact spaces attract 
much attention. On compact spaces the Gauss constraint restricts physical states into 
gauge invariant form, and due to this fact the theories are in a confinement phase at low 
temperature, and undergo a deconfinement transition at a critical temperature. Moreover, 
the large N gauge theories may have their dual description in terms of string theory on 
an asymptotic Anti-de Sitter (AdS) space pQ. For example, the partition function for 
a large N gauge theory on I x S 3 was computed in (2j [3], and it was shown that the 
partition function is of order 0(1) at low temperature and of order 0(N 2 ) above a critical 
temperature. In the dual gravity theory, the phase transition corresponds to the Hawking- 
Page transition [HE], where the thermal AdS space is dominant at low temperature and 
the AdS-Schwarzschild black hole is dominant at high temperature. 

In this paper, we study the thermodynamics of M = 4 U(N) super Yang-Mills theory 
We construct the orbifold theory in the following way. The manifold S 

1 The orbifold theory at zero temperature has been studied in [B], see also 0[H]. 
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has a U(l) symmetry along a cycle, and the orbifold is constructed by dividing rotation 
along the cycle. The background does not include any fixed point, but has a non-trivial 
cycle due to the orbifold procedure. We can introduce flux along the non-trivial cycle, 
which gives non-trivial holonomies to the fields. Therefore, the theory admits many vacua 
labelled by the choice of flux, and this makes the phase diagram richer. Along the non- 
trivial cycle, we can assign periodic and anti-periodic boundary conditions to fermions 
along the non-trivial cycle, and this leads to supersymmetric and non-supersymmetric 
theories at zero temperature. We study the gauge theories perturbatively with respect to 
the 't Hooft coupling A = iVgy m §| In this paper, we consider the zero 't Hooft coupling 
limit, where the theories reduce to free field theories. Even in this limit we observe phase 
transitions due to the compactness of the base manifold. 

In this orbifold case the dual gravity description is also available, thus the phase di- 
agram can be extended into the strong coupling region. For the case with the periodic 
boundary condition, the dual geometry is the orbifold of the thermal AdS space or its 
deformation by localized massless states in the low temperature phase. In the high tem- 
perature phase, the dual geometry is the orbifold of the AdS-Schwarzschild black hole or 
its deformation. If we perform the T-duality along the non-trivial cycle, then we obtain 
k NS5-brane configuration, which is parametrized by the positions of k NS5-branes [BJ. 
For the case with the anti-periodic boundary condition, there are localized tachyons at 
the fixed point of the thermal AdS orbifold in the low temperature phase. The conden- 
sation of localized tachyon may resolve the orbifold singularity like [9] and lead to the 
deformed geometry obtained in [TUl [TT] called as Eguchi-Hanson soliton. In other words, 
the gauge theory gives the dual picture of the localized tachyon condensation discussed 
in [9jlf| At enough high temperature, the dual geometry should be the orbifold of the 
AdS-Schwarzschild black hole, and there are no localized tachyons in the geometry. 

Following the analysis in [21 [3], we obtain the partition functions for the gauge theories 
in terms of a matrix integral. At low temperature, we find the leading contribution comes 
from the Casimir energy of the theories on S 3 /7*k- For the case with the periodic boundary 
condition, the Casimir energy is the same for all the vacua and the same as the mass of 
the thermal AdS orbifold. For the case with the anti-periodic boundary condition, the 
Casimir energy is smallest for the vacuum dual to the deformed geometry. Interestingly, 
the Casimir energy is roughly 4/3 times the mass of the Eguchi-Hanson soliton. In the 
high temperature limit, the partition function behaves in the same way for all possible 
holonomies and spin structures. Near the critical temperature, we can perform an analytic 
computation by using the Gross- Witten ansatz [T7j as an approximation, and we can 

2 On the compact space, we have a tunable dimensionless parameter RA, where R is the radius of S 3 
and A is a cut off scale. If we take RA <C 1 , then the Yang-Mills coupling can be set small even at low 
energy. Even with this fact we set R = 1 for simplicity. 

3 Previous attempts to apply the AdS/CFT correspondence to the localized tachyon condensation have 
been given in [TJ HI HI ED3 EES] • 
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discuss the dominant contribution to the total partition function. 

The organization of this paper is as follows. In the next section, we define the orbifold 
gauge theories on R x S 3 /7*k with generic holonomy by utilizing the standard orbifold 
method as in [18]. We compute the partition function for the orbifold gauge theories 
by following [21 E]- In section [3] we analyze the partition function and discuss the phase 
structure. We observe that the dominant contribution comes from the Casimir energy at 
low temperature, and we compute the Casimir energy for the case with general holonomy. 
Near the critical temperature we solve the partition function analytically by making use 
of the analysis in [T7j. At very high temperature, we obtain the partition function as an 
expansion of the temperature T, which does not depend on the choice of holonomy. In 
section H] we analyze the large 't Hooft coupling limit in the dual gravity description. We 
discuss the relations to the arrangement of k NS5-branes for the case with the periodic 
boundary condition and to the localized tachyon condensation for the case with the anti- 
periodic boundary condition. Section is devoted to conclusion and discussions. In 
appendix A we compute the partition function of single scalar particle as an example. 
In appendix B the index proposed in [19J is computed for the supersymmetric orbifold 



2 Orbifold gauge theories 

We consider M = 4 super Yang-Mills theories on R x S^/Z^ with large N U(N) 
gauge symmetry. The 3-sphere S 3 has £0(4) ~ SU(2)x x SU(2) 2 isometry, and we 
divide the gauge theories by 4n/k rotation along the x-cycle of U(l) x C 5 , f/(2) 2 . Since 
7Ti(S 3 /Zk) = Zfc, we can assign a non-trivial holonomy matrix V = Pexpi§A x along 
the non-trivial cycle. Using the U(N) gauge transformation, we can set V as a diagonal 
matrix V = diag(f2i, ■ • • , Ojv)- Because of the condition V k — 1, the element should be 
a k-th root of unity ttj = 1. Therefore, we can label the vacua by (no, • • • , nk-i) with 
=o n i, where nj represents the number of j such that Qj — u 1 (u — exp(2ni/k)). 
Then the orbifold theories are defined by projecting the Hilbert space into the orbifold 
invariant subspace as below. In the following we consider two special vacua. One is 
the Zfc symmetric vacuu with m = N/k for all /. Due to the symmetry the dual 
geometry can be identified as the standard orbifold. The other vacua do not preserve the 
Zfc symmetry, thus the dual geometry should be a deformation of the orbifold. The other 
important vacuum is with the trivial holonomy V — 1 or equivalently n Q = N. In this 
case the Z^ symmetry is maximally broken. 

4 Similar computations were done in [2U] for quiver gauge theories, which are constructed as orbifolds 
different from ours. 

5 In this case N is assumed to be TV = kZ. However, for large N and finite k, the difference from the 
general N case should be negligible. 
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2.1 Partition function 



We would like to compute the partition function of gauge invariant operator in the 
orbifold theories. In general, the counting of gauge invariant operator is an involved task. 
Fortunately, it was shown in [21 E] that the partition function of gauge invariant operator 
can be written in terms of single-particle partition functions as 



Z(x)= J [dU] exp 



OO 1 

. K n=l n 



(2.1) 



At this stage, a gauge group G could be arbitrary, and the sum is taken over all represen- 
tation 1Z of the gauge group G. We denote Xn(U) as the character for representation 7Z, 
and [dU] as the Haar measure for the group element U. The partition function of single 
particle in the representation TZ is computed by 

z*{x) = (2-2) 

E 

where E denotes the energy eigenvalue. The condition of gauge invariance comes from 
the integral over U, and the variable U will be identified as the holonomy matrix along 
the thermal cycle. 

We can adopt any gauge group G and representation TZ in the formula (12.1 1) , and several 
interesting examples have their dual gravity description. The most famous one arises from 
the iV D3-brane worldvolume theory, which is dual to superstrings on AdS§ x S 5 . The 
theory is M = 4 super Yang-Mills theory onMxS 3 with gauge group G = U(N), where the 
states are in the adjoint representation. Our interest is on the orbifold gauge theories with 
holonomy along the non-trivial cycle, where the existence of the holonomy (no, • • ■ , nk-i) 
breaks the gauge symmetry into G = ]l/=o^(^/)@ With respect to the broken gauge 
group, the states are in the adjoint representation for U(ni) or in the bi-fundamental 
representation (nj,nj) for U{nj) x U(nj). 

The spectrum of the orbifold theories can be obtained by projecting the spectrum 
on S 3 into the orbifold invariant subspace. The spectrum on S 3 can be obtained from 
the spherical harmonic analysis as in |21j. The theory includes 6 scalers, a gauge field 
and 4 Majorana fermions. The scalars can be expanded by the scalar spherical harmon- 
ics Sj jnij fh(Q), where Q represents the coordinates of S 3 . The eigenfunctions of Laplace 
operator on S 3 are given as 

V 2 ^,™,^) = -](] + 2)S J<m! rn(n) . (2.3) 

The labels (m,m) are eigenvalues of J 3 and J 3 for SU(2)i and SU(2)2, and they run 
— j/2, — j/2 + !,•■■ ,j/2 — l,j/2. The projection into the orbifold invariant modes is 



6 The diagonal U(l) parts of each U(ni) may be decoupled from the rest, but the difference can be 
ignored when ni are very large. See, however, appendix B for the case of an index. 
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performed by the projection operator P = ~ Y^iZq T 7 , where T represents the orbifold 
action. For a bi-fundamental state (rij, rij) or an adjoint state with I = J, the orbifold 
action is given by 

where u = exp(2iri / k) . The orbifold action consists of two parts. The first part is the 
phase shift due to the rotation along the x-cycle. The second part is the holonomy 
for the bi-fundamental state (nj,nj). We can see from (12 .4p that the orbifold invariant 
modes are restricted to 2m = J — I mod k. As a notation we define < L < k subject 
to L = J — I mod k. Now we can compute the partition function for the single scalar 
particle (12. 2p as 



z s ( x ) 



+ 



(x L+1 - 


- x L+3 + X 






(1 -X 2 ) 2 ! 


'l-x k ) 2 


(L + l 


)x L+1 - (L 


- l)x L+3 — (L — l)x k - L+1 + (L + l)x k - L+s 


(1 -x 2 ) 2 (l -x k ) 



(2.5) 



See appendix A for the detail. We have used the fact that the energy is given by E = j + 1 
for a scaler on S 3 conformally coupled to gravity. 

We move to the gauge field, which is expanded by the vector spherical harmonics 
^m,m(^)0 We use the notation such that the vector index is contracted with an auxiliary 
unit vector £ M as K^j^.„£ M . The vector spherical harmonics V^ mfh and V~ mfh belong to 

2 ' 2 

Laplace operator on S 3 are 



the representations (31,32) = (4^i V) an< ^ ("Stj 2 ^ - )) respectively. The eigenvalues of 



V 2 V^(ft) = -(j + l) 2 V± m<fh (n) . (2.6) 

The orbifold action to the vector spherical harmonics is the same as in the scalar case 
(12 .4p . since the vector index is contracted with an auxiliary unit vector. Therefore, the 
orbifold projection allows only the modes with 2m = J — I mod k for the bi-fundamental 
state with (ni,nj). The partition function is then given by 

ZV+[X) = (1-Z 2 ) 2 (1-X*) 2 

(2.7) 

(L + 3)x L+2 - (L + l)x L+A - (L - 3)x k - L+2 + (L - l)x k ~ L+4 
+ (1 -x 2 ) 2 (l -x k ) 



A longitudinal mode is expanded by the scalar spherical harmonics as VS 1 , and we do not consider 

it. 
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f o r V 3 \m,fh and 



IyJ (x L - x L+2 + x~ L - x~ L+2 )kx k 

Z v~^ x) = (1 -x 2 ) 2 (l-x fc ) 2 

(2i 

(L - l)x L — (L — 3)a; L+2 - (L + l)x k ~ L + (L + 3)x k ~ L+2 



+ 



x 2 ) 2 (l — x k ) 



for Vf mm . We have defined L = J — I mod (0 < L < k) as before. 

Fermions are expanded by the spherical harmonics F^ mfh (fl) and Fj~ mm ($V), which 
belong to (j'1,.7'2) = (|> ^f - ) an< ^ (^i - ;!) - ^ ne s Pinor index is again contracted with an 
auxiliary spinor x° as ^>i,m;aX a ' The eigenvalues of Laplace operator on S 3 are 

V 2 i^(fi) = -(j + \) 2 F± m ^) . (2.9) 

For fermions we can assign two types of boundary conditions along the non-trivial cycle. 
The orbifold action depends on the boundary condition as 

T = ± e 4mj3/k u I - J , (2.10) 

where + and — means the periodic and the anti-periodic boundary conditions, respec- 
tively. For the periodic boundary condition, the orbifold invariant modes are given by 
those with 2m = J — I mod k. The anti-periodic boundary condition can be assigned 
only for even k, and the restriction is shifted by k/2 as 2m = J — / + k/2 mod k. For the 
periodic boundary condition, the partition function can be computed as 

i j (x L+ 5 — x L+ i + — x~ L+ ^)kx k 

Z F+\ X ) 



' ' 1 " (1 ~.r 2 )-'(l - r h )- 

(L + 2)x L+ l - Lx L+? 2 — (L — 2)x k - L+ l + Lx k - L+ i 



[2.11) 



+ 



(1 -x 2 ) 2 (l -x k ^ 



for F j,m,rh and 



Z TP — I X 



x L+ 2 — x L+ 2 4- x L+ 2 — x L+ 2)kx 



I ■ ( 1 ) 2 (1 x k ) 2 

Lx l+\ _ ( L _ 2 )x L+ i - Lx k - L+1 * + (L + 2)x k - L+ i 



:2.i2i 



{1 - x 2 ) 2 {l - x k ) 

for F~ m fh . For the anti-periodic boundary condition, we should use L = J — I + k/2 mod 
k with < L < k instead of L = J — I mod k. 

Now we can write up explicitly the partition function (12.11) for M = 4 super Yang- 
Mills theories onlx 5 3 /Zfc with holonomy (no, • ■ • , rik-i)- The total partition function 
is given by summing over all the vacua. It is useful to use the formula for the character 
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of bi-fundamental representation as X(ni,nj){U) = Tr Uj Tr Lf\, where the trace is taken in 
the fundamental representation. The partition function is then given by 



Z(x) = J [JJdE/j]exp 



OO 1 

EE-4' J (^)Tr(f/7)Tr(^) 
tj n =l n 



(2.13) 



1 

where the single-particle partition function is summarized as 

Z y(x) = 64< J (x") + zftix") + 4 J -{x n ) + (-l) n+1 4[4((x n ) + z'/ix")} (2.14) 

for the case with the periodic boundary condition and 

JJ(x) = Qz^ J (x n ) + z\}{[x n ) + z^t{x n ) + (-l)" +1 4[4+ +f (x n ) + z F J - + ~ 2 (x n )} (2.15) 

for the case with the anti-periodic boundary condition. 

Finally let us remark on the difference from the D-brane worldvolume theories localized 
at the fixed point of C n /Y with n = 2,3 [IE]. Since the orbifold action acts trivially to 
the worldvolume in those cases, only bi-fundamental matters with 1,1 ± 1 and adjoint 
gauge fields (and matters) are left under the orbifold projection. On the other hand, the 
orbifold action rotates S 3 by 4ir/k in our case, there is a m dependent phase in (12.41) . 
Due to this effect, bi-fundamental states with every pairs of /, J (and adjoint states with 
I = J) survive the projection each for matters, gauge field and fermions. The difference 
would be significant if we compare our case with the duality between superstrings on 
AdSc, x S 5 /Y and the gauge theory coming from D3-branes at the fixed point of orbifold 
action Y [22]. 



2.2 Path integral formulation 

In the previous subsection, we have obtained the partition function of gauge invari- 
ant operator (12.131) in terms of integral over the group manifolds. However, we cannot 
determine the overall pre-factor in the formulation. In this subsection, we re-derive the 
partition function in the path integral formulation. In this derivation, we obtain the nor- 
malization depending on the Casimir energy of the gauge theories on R x S 3 fZk- The 
Casimir energy will be important when we consider the phase structure at low tempera- 
ture. Moreover, we can identify Uj as the holonomy matrix for U(ni) gauge group along 
the thermal cycle. 

The path integral for the partition function with a finite temperature T may be com- 
puted on S 1 x S 3 /Zk, where S 1 is the thermal cycle with periodicity (3 = 1/T. Along 
the thermal cycle we assign the anti-periodic boundary condition for the fermions§| We 

8 Due to this boundary condition, supersymmetry is always broken at a finite temperature even for 
the theory supersymmetric at zero temperature. 
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start from fixing the gauge symmetry and then introduce the Faddeev-Popov determinant 
conjugate to the gauge fixing. We adopt the Coulomb gauge 



V a A a = (2.16) 

with Vffl as covariant derivatives along the S 3 direction (a = 1, 2, 3). If we do not include 
a non-trivial holonomy, then there are spatially constant modes of the gauge field. The 
presence of holonomy (n , • ■ ■ , n k -\) breaks the gauge group into Hi U(nj) and spatially 
constant modes are left only for Il/^( n /)J^I The time-dependence of these modes is not 
fixed by the Coulomb gauge (12.161) . and we fix these degrees by 



d ' a ' = ' « '= voWm J daA '" (2 ' 17) 

where the integration is performed over S^/Z*;. 

First we consider the Faddeev-Popov determinant conjugate to (12.171) . which is given 

by 

A FP = det'idtDl) , D\ = d t - {[a 1 , } . (2.18) 

The determinant is taken over the non-zero modes. Diagonalizing the zero modes as 
a 1 = diag(af , • ■ ■ a£ ), the measure can be written as 

da 1 = Y[daj J] \a{ - a]\ , (2.19) 

where the Van der Monde determinant arises from the integration over the off diagonal 
elements. Now that the bosonic modes are periodic along the thermal cycle, they can 
be expanded by the function exp(2irint / (3) with n £ Thus the determinant can be 
written as 



. t T-r T-r 27rm / 2irin , , \ 

a fp = n n — r — r - ^ - a ^ ) • ( 2 - 2 °) 

i,j n^O P V P 



With the help of the formula E[^=i(l — x 2 /n 2 ) = sin7rx/(7rx), we find up to an overall 
factor 

[dUj] = da'A 1 ^ = TJ dai J] sin 2 f ^~°^A ; ( 2 .21) 

i i<j \ ^ J 

which is the Haar measure of Ui = exp(i/3o/). 



9 The projection under the orbifold action (|2.4|) removes spatially constant modes with m — for 
I ^ J sectors. 

10 For fermionic modes, we should replace n by n + 1/2 since we have assigned anti-periodic boundary 
condition along the thermal cycle. 
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The Faddeev-Popov determinant conjugate to (12.161) is given by 

detV a £> a = J[dcdc]exp(-cV a D a c) , (2.22) 

which should be added to the action of M = 4 super Yang-Mills theory. Notice that 
the ghosts are expanded by the scalar spherical harmonics projected into the orbifold 
invariant modes. After integrating over the massive modes including the c-ghosts, the 
partition function is given in terms of integral over the zero modes as 

Z ( T ) = !\R dU Ae~ S{U) ■ (2-23) 
J i 

Let us first compute the contribution from the gauge field and the c-ghosts. Since the 
longitudinal modes V S and A t (except the zero modes a 1 ) are expanded by the scalar 
spherical harmonics, the contributions to the path integral from the c-ghosts and the lon- 
gitudinal modes cancel out. Therefore, the contribution reduces to the Gaussian integral 
over the vector spherical harmonics, which is evaluated as 

S(U) = WY.{4i(E)+n I v J (E)}\ndet(-D? + E 2 ) . (2.24) 

Z I,J E 

We denote ny± (E) as the degeneracy of eigenstates with E in the representation (rij, fij). 
Following the computation in [3], we find 



S(U) = Wf3n I n J Y.[n v {{E)+n v i{E)}E 

/ I,J E 



(2.25) 



OO 1 

E E 34 J + (e- n/T ) + 4-(e- n/T )] T*W) T±(U] 



I.J n =l 71 



In the same way, we can compute the contributions from scalars and fermions and sum- 
marize all the contributions aa^ 

OO 1 

S(U) = /3V - E E -4-V 1/T ) Tr(C/?) MU\ n ) , (2.26) 

I, J n=l n 

where z^ J are the single-particle partition functions (12.141) or (12.151) . The first term is the 
Casimir energy 

V o = \Y. n i n jY.\Sn/{E)+n^i{E)+n I ^{E)-An 1 ^ . (2.27) 

1 I.J E 



1:L In the expression of (|2.23[) with (|2.26p . the normalization has been set by dividing holonomy in- 
dependent factors. Introducing the holonomy (no,-- - ,nj) along the x-cycle, the lowest modes in 
U(N)/ Yii U(ni) become space-dependent, and hence they can be fixed by the Coulomb gauge (12.161) 
instead of (|2.17[) . The normalization of each Faddeev-Popov determinants may depend on the choice of 
holonomy along the ^-cycle, but the sum of both should not. 
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Compared with the expression of (12.131) . the integral variables Uj are identified with the 
holonomy matrices Ui = exp^ifia 1 ) with respect to the gauge group 11/ U (nj). In this way 
we can see that the previous expression only includes the finite temperature contribution. 
The zero temperature contribution, which comes from the Casimir energy, should be 
included in the partition function. 



3 Phase transitions of the gauge theories 

In the previous section, we have obtained the partition function of gauge invariant 
operator in terms of integral over Uj as (I2.23H with (12. 26ft . In this section, we perform 
the Ui integral in the large N limit] 12 ! and examine the phase structure of the orbifold 
theories. For large N and fixed k, it is natural to assume that rij in the label of holonomy 
(no, • • • , Jifc_i) are very large. In case that some of nj are very small, then they may 
be set zero in this limit. We consider two specific vacua in the following. One is the 
Zfc symmetric holonomy vacuum with nj = N/k and the other is the trivial holonomy 
vacuum with n = N. In these cases our assumption is valid. In the next subsection we 
study the low temperature phase, and in subsection 13.21 we focus on the Casimir energy 
contribution. In subsection 13.31 we obtain an analytic expression under an approximation 
near the critical temperature. In subsection 13 .41 we take the high temperature limit, where 
the analysis becomes simpler. 



3.1 Critical temperatures 



It is convenient to diagonalize the eigenvalues of holonomy matrix Ui as exp(i8i^) 
with -7i < Qn < 7T@ For large nj the discrete elements may be replaced by a continuous 
parameter Oj with a density p T (9i). The density has to satisfy p J (0i) > and the nor- 
malization is set as J^ n p 1 (9j)d9j = 1. In this approximation the effective action (12.261) 
becomes 



51/(0/)] = /3V -Y^^inj / d9 I d9' J p I (9 I )p J (9' J ) 
i, j J 



1 



Srj In 



6', 



sin 



(3.1) 



+ £-4' J (*)cos(n(0 J -^)) 



n=l 



n 



with x = e 1//T . The first term in the bracket arises from the change of measure [dUi] — > 
[d9i t i\. In terms of the Fourier transform p I n = J d9ip I {9i) cos(n#/)Jl3 the effective action 



12 Finite N effects may be examined by following the analysis in [551 [Ml [25] 
13 The eigenvalue Oi^ is related with the zero modes in (|2.17[) as Oj : i = (3a{. 
14 We assume that 8] is distributed symmetrically around 6j = 0. 
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( 13. lp is given by 



where 



S[ P i] = PV + Y,nmj £ pipiV^ \x) , (3.2) 

I,J n=l 



Vy(x) = -{Si,J - z 'n J ^)) ■ (3-3) 
n 

At enough low temperature, the repulsive force coming from the first term of (13.31) 
dominates, and the uniform distribution = for n > 1 is the classical solution to the 
effective action (13. 2p . Therefore, there are no order 0(N 2 ) (nor order O(N)) contribu- 
tions from the effective action except for the Casimir energy term /3Vq. An order 0(1) 
contribution comes from the Gaussian integral as 

OO ^ 

n 777 7W7Y ' (3 ' 4) 

n =i det [n^jVn (x)J 

where the determinant is over the labels < I , J < k. As the temperature increased, the 
second term of (I3.3P contributes to the potential, and the determinant would vanish at a 
critical temperature x c = exp(— 1/T C ). Above the critical temperature, the distribution 
becomes non-uniform and the classical contribution is of order 0(N 2 ). 

Let us examine two concrete examples. We start from the trivial holonomy case, where 
the action (13. 2p reads 

OO 1 

%°] = W + N 2 J2 IpT-(1 - #°(*)) • (3-5) 

n=l U 

At enough low temperature, the coefficients of |p°| 2 are positive, and p° = for n > 1 is 
the saddle point. Now that the coefficients are lxl matrices, the determinant (13.41) is 
simply 

1 



n=l 1 



Z°n Q (x) 



(3.6) 



We have changed the normalization such that only the Casimir energy term is left at 
x = 0. As the temperature increased, the coefficients of |p°| 2 become smaller, and at 
a critical temperature T = T c , a coefficient vanishes. Since the single-particle partition 
function is a monotonically increasing function of x, the first zero comes from the n = 1 
part when 1 — Zi'°(x) = 0. The critical temperatures x c and T c are summarized for small 
k in Table [U We should note that for k = 1 the critical temperature reduces to the one 
for Rx S 3 case E] as x c = 7 - 4^ = 0.071797 or T c = 0.379663. 

Another interesting case may be with the symmetric holonomy rij = N/ k for all /. 
In this case the action ( 13.21) is given by 

N 2 °° 1 

%n] -W + ¥ EE rtPi-Vv - 4' J (*)) • (3.7) 

I,J n=l 71 
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k 


x c (periodic) 


T c (periodic) 


x c (anti-periodic) 


T c (anti-periodic) 


2 


0.095663 


0.426090 


0.095663 


0.426090 


4 


0.104448 


0.442661 


0.127999 


0.486445 


6 


0.104684 


0.443104 


0.139545 


0.507777 


8 


0.104689 


0.443113 


0.142528 


0.513290 


10 


0.104689 


0.443113 


0.143136 


0.514414 



Table 1: The critical temperatures x c = exp(— 1/T C ) and T c for the Af = 4 super Yang- 
Mills theories on M x S 3 /^ with the trivial holonomy V — 1. We set k = 2,4,6,8, 10. 
The periodic and anti-periodic boundary conditions are assigned for the fermions along 
the x-cycle. 



Notice that the coefficients of p ! n take the form of a circulant matrix since z^ J only 
depends on the difference I — J. Using the formula for a circulant determinant (IB. 111) , 
the determinant ( 13. 4ft can be written in a compact form as 

oo k— 1 -y 

II II ; u o,j, ; • ( 3 - 8 ) 

n=l 7=0 1 — 2^ J=o w 2 « \ x ) 

As before it is enough to focus on the n = 1 factors. If we increase the temperature, then 
the denominator diverges when J2j z i' J { x ) = 1- Among the other factors, this factor gives 
the divergence with smallest x c since z®' (x) is positive for all J. Remarkably the critical 
temperature is the same for all k and for both the spin structures as x c = 0.071797 or 
T c = 0.379663 as in the K x S 3 case. Actually this is an expected result since the sum of 
all sectors with weight one J2j z i' J ( x ) is the same as the single-particle partition function 
for the M x S 3 case. 



3.2 Casimir energies 

At low temperature, the determinant ( 13. 4ft is of order 0(1) and the contribution from 
the Casimir energy (12.27P is dominant. The Casimir energy is an important quantity 
since it is supposed to correspond to the mass of the dual geometry. In order to compute 
Casimir energy, we have to sum over infinitely many states, and this may lead to a 
divergent result. Thus we have to choose a regularization, but it is a subtle problem 
for quantum field theory on a curve background [26]. Fortunately, we will find that the 
Casimir energies in our case are finite, thus we do not need to worry about this difficult 
issue. 

In order to regularize the infinite sum in the Casimir energy ( 12. 2 7ft . we first introduce 
a cut off factor e~ E ^ as 

£(6n 7 / + n# + rfyL - 4n# - An^Ee-^ , (3.9) 

E 
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and finally take the limit e = — > 0. This regularization may be justified by the fact 
that no divergent terms are left in the final results. The above regularization can be 
realized by using the identity (A = S, V ± ) F ± ) 

£^ J (£)£e-^ = -^4V e ), (3-10) 

E a€ 

so we need the expansion of single-particle partition functions by e up to 0(e) as 

rj 2 7e ek ek 3 ehik - L) eL 2 (k - L) 2 
Zs (e > ~ i*jfc ~ 180* + 36 + 360 6* 12*~ ' (3 ' U) 



u, - £ s 2 1*1 L 2 2e 5e* ek 3 

^' ±(e ) ~^k-7k ± 6 ± Tk TL± T + {—) 5L > + 45k-^6 + 360 

, 5eL(k-L) eL 2 (k-L) 2 

+ 6k 12k ' {6 - U) 



2 1 * 1 L L 2 83e e* e* 3 



e 3 * 4e* 12 1 12* 2 2* 2880* 72 360 
eL(k-L) eL 2 (k-L) 2 

+ ^I2F J laH-' (3 ' 13) 

Below we discuss the cases with the periodic and the anti-periodic boundary conditions 
separately. 

For the case with the periodic boundary condition, we have 

]T(64' J + rfy + n y_ - fag _ 4r^)£ e -^ = A + ( e ) . ( 3 .14) 

As mentioned before, the final expression does not depend on the cut off parameter e, 
which might be due to the large supersymmetry. Moreover, the Casimir energy does 
not depend on the indices J, J, and this means that the Casimir energy is the same for 
all choices of holonomy (no, • • • , n^-i). This is consistent with the argument in [6] that 
the vacua with different holonomy are degenerated at zero temperature. Notice that the 
Casimir energy 

V = JVA (3.15) 

is precisely the same as the mass of AdS^fLk as we will discuss below. 

For the case with the anti-periodic boundary condition, the above cancellation among 
/, J dependent terms does not occur in general. With L = J — I mod * (0 < L < k) we 
find 

£(64< J + n T / + + rfyL - 4n F J + H - An F i +i )Ee^ E 

(3.16) 



3 



3 * * 
8* + 6 "24 +L 



* - - ( 2L + - 
3 V L 



+ 0(e) 
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for < L < | and L is replaced by k — L for | < L < k. Notice that the divergent 
terms proportional to 1/e 4 and 1/e 2 cancel out even in this case. The Casimir energy 
depends on the choice of holonomy due to the L-dependence of the zero point energy 
(I3.16p . The dominant contribution to the total partition function comes from the vacuum 
with smallest zero point energy, which is realized with the trivial holonomy. This is 
because the value inside the bracket in (I3.16P is always positive for all < L < | if we 
set k > 4@ The Casimir energy for the trivial holonomy case is given by 

^^(s4-s)' < 3 - 17 > 

which will be compared with the mass of the dual geometry (14.121) . Another interesting 
case may be with the symmetric holonomy. In this case, we sum up every L with the 
same weight, and hence we except that the cancellation between the sectors with L and 
L + k/2 occurs. This can be confirmed by a direct computation, and the Casimir energy 
is obtained as fl3~T51) l"j 

Before finishing the arguments on Casimir energy, we would like to make a comment 
on the validity of regularization adopted here, even though the divergent terms cancel 
out in the final expressions. Let us write the radius R of S 3 explicitly such as E = 
(j + 1)/R for the scalar case. Then the divergent parts are proportional to p 4 R 3 and 
p, 2 R in each single-particle partition function. For quantum field theory on a curved 
background, the divergent terms of energy momentum tensor should be absorbed by the 
renormalization of coefficients in the Einstein-Hilbert action. In our case the divergent 
terms may be absorbed by the counter terms a/i 4 / y/g and bp 2 j ^fgTt. See |26j for more 
detailed discussions. 



3.3 Just above the critical temperatures 

The eigenvalues distribute uniformly due to the repulsive potential at low temperature, 
however the eigenvalues get together due to the attractive potential above the critical 
temperature. In particular, the densities may be gaped and vanish except for — 9j c < 
0i < @ic- The condition that an eigenvalue 9[ does not feel any force is obtained as 

„ /a ft' \ ^— 1 oo 

n] J dO'jp^e'j) cot ±—L = 2 £ njnj £ z l n \x)p{ sin(r^) (3.18) 

\ / J=0 n=l 

from the action (13. ip . The general solutions subject to the normalization condition p^ = 1 
can be obtained by following [271 E] in principle. However, the analysis is quite compli- 
cated generically, so we adopt an approximation by setting z^ J (x) = for n > 1. This 

15 For k = 2 the cases with the periodic and anti-periodic boundary conditions lead to the identical 
result. 

16 The Casimir energy in this case was already computed in [101 111] by following the general method of 
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approximation may be justified for small x ~ x c as in Table [T] by the fact that (x) 
with n > 1 is much smaller than z[ ,J {x). In the following we will explicitly solve these 
equations for the trivial holonomy case with V = 1 and the symmetric case with 
nj = N/k. 

Let us begin with the trivial holonomy case, where the condition (I3.18P reduces to 

J d6' p°(6' ) cot (^^j = 2z 1 '°(x)p° l sm6 . (3.19) 

This case is almost the same as the IxS 3 case analyzed in [21 [3] . The solution is given 
by the form of the Gross- Witten ansatz jTTj as 



P (0o) = . 2e c \ sm 17 - sm 77 cos — (3.20) 
7r sin I y 2 2 2 

for — 6 C < 6q < 9 C and zero for otherwise. The parameter 9 C satisfies 



6c 
2 



■ 2 u c 1 

sin — = 1 



(3.21) 



With this solution we can compute the classical action and the free energy F = —T In Z 
T(S) as 



w^ v °- T bk + 5" lsm2 |-il • < 3 - 22 > 



2 



Near the critical temperature, it is given as 

■j^ — Vq — -£-(T — T c )±z°Ae-V T )\T=T c + 0((T- T c f) . (3.23) 



For the purpose of comparison with the symmetric case, we draw plots of the free 
energies for k = 4, 6 in Figure [TJ In the Figure we have shifted the zero point energy by 
3/(16Jfe). 

For the symmetric case, the condition (13. 18ft becomes 



J de'tp 1 ^) cot ( 6 - L ^ L ) =2j24> J (x) P J n sme I , (3.24) 



and the generic solutions are quite complicated. However we only need the solution 
responsible to the phase transition at the critical temperature x c satisfying Sj^' J (x c ) = 
1. With the help of the symmetry, we assign that the densities of eigenvalue take the 
same form as p 1 ^ 1 ) = piO 1 ) for all /. Then the solution can be easily found as 



\ I n9 c Q 6 J 6 

piO,) = T— f| f in - - sin - cos - (3.25) 
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F/N"2 



periodic k-4 



. JTS~\ .5 0.55 0.6 0.65 



-1 L 
F/N"2 



-0.8 



periodic k= 



t4&\ 0.5 0.55 0.6 0.65 



F/N"2 anti-periodic k=4 

1 r 



).45 0~.~5 9^5-5- 



F/N"2 anti-periodic k= 

1 r 



).45 0.5 ~-~0t55 0-£._ 0.65 



Figure 1: Free energies F(T)/N 2 as functions ofT in the cases with the periodic and the 
anti-periodic boundary conditions and with k = 4,6. The solid lines are for the trivial 
holonomy case and the dotted lines are for the symmetric case. 



for — 9 C < 9 < 8 C and zero for otherwise. In this case 8 C satisfies 



• 2 1 
sin — = 1 — 



(3.26) 



This is valid for x satisfying J2j z i' J ( x ) — 1 ; and the equality holds for the critical tem- 
perature x = x c . The free energy is 



N 2 



Vn- 



t f i i . 2 e c i 

— 5-3- H — In sm 

fc \2sin 2 % 2 2 2 



(3.27) 



and near the critical temperature 



N 2 



T=T C 



+ (T-T c 



(3.28) 



See Figures [T] for k — 4,6. 

Let us compare the free energies for the above two cases. For the case with the 
periodic boundary condition, the free energy for the symmetric case is lower near the 
critical temperature, since the critical temperature is smaller in this case. Near the critical 



temperature, the free energy is proportional to T — T c as in ( 13.231) and ( 13. 28ft . and the 
coefficients in (13.281) behaves like i/kE From this reason the free energy for the trivial 
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This is checked for small k by a numerical computation. 
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holonomy case is lower at slightly higher temperature for large k as in Figure [TJ For the 
case with the anti-periodic boundary condition, the free energy for the trivial holonomy 
case is lower than for the symmetric case due to the Casimir energy Vq. Since the 
Casimir energy behaves like k 3 , the difference of free energy becomes bigger for larger k. 



3.4 High temperature behaviors 



At higher temperature the approximation in the previous subsection is not valid any 
more, and the contribution from n > 1 terms in (13.181) should be taken into account. 
Fortunately, the analysis becomes simpler when we take the high temperature limit T ^> 1. 
This limit is the same as the limit of large radius R of S 3 , where we can use the flat 
space approximation. In this limit, the densities of eigenvalue may be given by the delta 
function, thus we set p 1 ^ = 1 for all n. From f 1 3 . 3 H we find 



E Vi = -C(4)T 3 



n=l 



16 
T 



+ 1 



1 



16 

T 



C(2)T 



2 
k 



O 



(3.29) 



which may be read off from the expansion of z£ (e _e ) by e. The coefficients of T 3 terms 
can be given by the degrees of freedom in the flat space limit, where the volume should be 
divided by k. Since the above expression does not depend on L, the free energy behaves 
in the same way for all the vacua and for both the spin structures in the high temperature 
expansion. This is quite natural since the high energy excitations should not depend on 
the vacuum structure. In summary, the free energy F = —T In Z and the expectation 
value of energy E = — In Z are given by 




7T 



T 2 + 0{l) 



E 




71 1 



7T 



T 2 +0(1) (3.30) 



for every choice of holonomy and for both the spin structures. This energy was already 
computed in [28] in a different way. 



4 Dual gravity description 

In the limit of large 't Hooft coupling the dual gravity description is more appropriate 
to discuss the phase structure of the gauge theories. In the dual picture the confine- 
ment /deconfinement phase transition is described by the Hawking-Page transition [HE]. 
In this section we investigate the Hawking-Page transition in the dual geometries and try 
to see whether the phase structure continues from the one at zero 't Hooft coupling. In the 
next subsection we review the Hawking-Page transition between the thermal AdS space 
and the AdS-Schwarzschild black hole. In subsection 14.21 we move to our orbifold cases. 
With the periodic boundary condition, the geometries are given by various arrangements 
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of k NS5-branes in a T-dual picture. With the anti-periodic boundary condition, there 
are localized tachyons at the fixed point of the thermal AdS orbifold. The condensation 
of the localized tachyon is discussed in subsection 14.31 



4.1 Hawking-Page transition 

The boundary gauge theory at finite temperature may be defined on S 1 x S 3 with a 
thermal cycle, and we have to include all geometries with the same boundary condition to 
compute a path integral in gravity theory. The geometries are obtained by extending the 
boundary S 1 x S 3 into the bulk, which can be done in two ways. One of the geometry has 
the topology of S 1 x B 4 , where the 4 dimensional ball B 4 is obtained by filling the inside 
of S 3 . This geometry is the thermal AdS space, where the thermal cycle has a periodicity. 
The other geometry has the topology B 2 x S 3 , which is the AdS-Schwarzschild black 
hole. We will see that the thermal AdS space is dominant at low temperature and the 
AdS-Schwarzschild black hole is dominant at high temperature. 

The metric of the thermal AdS*, is given by 

dv 2 v 2 
ds 2 = g{r)dt 2 + — + r 2 ds% , g{r) = - + 1 , (4.1) 

where ds 2 s3 is the metric of S 3 and the Euclidean time is periodic t ~ t + The mass of 
the thermal AdS$ was computed in [29] as 

3tt/ 2 , s 

M = —-- 4.2 
32G 5 v ; 

by utilizing the boundary stress tensor method. This is precisely the same as the Casimir 
energy of M = 4 U(N) super Yang-Mills theory on R x S 3 given in (13.151) with k — 1. 
The 5 dimensional Newton constant G5 is written as G5 = ttI 3 / (2N) in terms of the dual 
gauge theory. For simplicity we set the AdS radius as I = 1. 
The metric of the AdS-Schwarzschild black hole is 

dv 2 v 2 
ds 2 = h(r)dt 2 + TTT + r 2 ds 2 s3 , h(r) = r 2 + 1 - 4 , (4.3) 

h[r) r l 

where the period of the Euclidean time is given by (3 h = 2nr + /(2r 2 + 1). We denote 
r + as the horizon satisfying h(r+) = or equivalently r 2 = r\ + r^cj The Mass of the 
AdS-Schwarzschild black hole is 



M = + — — , (4.4) 



18 There are two solutions r± to this equation, which implies that there are two types of black holes. 
The bigger and smaller ones r + ,r_ are the radii of horizons of big and small black holes, respectively. 
We only consider the big black hole since the small black hole has a negative specific heat and hence it 
is unstable. However the unstable saddle point may be important to understand the phase structure as 
pointed out in [25] . 
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where the second term is the constant AdS§ contribution. This mass can be expanded by 
T = l/(3h at high temperature as [28] 

3"7T 

M = ^r(n 4 T 4 - tt 2 T 2 ) + 0(1) , (4.5) 
8G 5 

which is about 3/4 times the energy given in (13.301) with k = 1q 

Now we can compute the partition function in the gravity description. In the classical 
approximation, the partition function is obtained from the classical actions for the geome- 
tries, which are proportional to the volume as S = V/{2ttG^). Since the volume diverges 
for both cases, we introduce a cut off r m and examine the difference. The Euclidean time 
periods are set equal at r m as fliJg(r m ) = f3h\Jh{r m ). Then the difference of the classical 
action is computed as [H [S] 

, V B H{r m ) ~ VTAds(r m ) n 2 (rl - rjj 

2ttG 5 4G 5 (2r 2 + 1) ' [ ^ } 



For small r + (low temperature) the above quantity is positive, which means that the 
thermal AdS space dominates the partition function. The phase transition occurs at 
r + = 1 or T = 3/(2-7r), and above the critical temperature the AdS-Schwarzschild black 
hole dominates. 

In the dual gauge theory, the Polyakov loop (Tr U) = (Tr P exp % § A t ) is an important 
order parameter. In the confinement phase, the Zjv symmetry of the theory is unbroken, 
and hence the Polyakov loop vanishes, which is realized by uniformly distributed eigen- 
values 0j. In the deconfinement phase, since the Zjy symmetry is broken, the Polyakov 
loop may have a non-trivial value and the eigenvalues are distributed non- triviallyH 
The Polyakov (Wilson) loop along the path C may be computed in the gravity side as 
exp(— A), where A represents the minimum area of the worldsheet with boundary C sub- 
ject to a regularization [2U E2]- For the thermal AdS space, the thermal cycle is not 
contractible, so the area is infinite and this leads to the vanishing Polyakov loop. For the 
AdS-Schwarzschild black hole, the thermal cycle shrinks at the horizon, so the area could 
be finite and hence the Polyakov loop can take a non-zero value. See [5] for more detail. 



4.2 Phase transitions of the orbifolds 

Due to the Z& symmetry the vacuum with the holonomy rii = N/ k for all / is supposed 
to be dual to the standard orbifold. At low temperature, the dual geometry is given by 
the orbifold of the thermal AdS space, whose metric is fl4.ll) with ds 2 s3 replaced by 

ds 2 OTh = - [{dx + cos dd<pf + dd 2 + sin 2 dd<p 2 } , (4.7) 

19 It was argued in [3H] that the origin of the 3/4 difference is the same as the one of [30], where the 
entropy of black 3-branes is compared with the state counting on D3-brane. 

20 Because the eigenvalues collapse in the Z^r symmetric way, the Polyakov loop may vanish after taking 
the average. We may use (| TiU\) or (TrC/ 2 ) to avoid this subtlety. 
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where the variables run < # < 7r, < <p < 2tt and < x < Air/k due to the orbifold 
identification. This metric of S l3 /Z/ C may be useful since we can easily see that the orbifold 
action acts on the x-cycle as x — > x + ^jk and that there is no fixed point on this space. 
The mass of the thermal AdS orbifold is 

M = , (4.8) 

which is 1/k times the mass of the thermal AdS space (14.21) . This should be compared 
with the Casimir energy with the symmetric holonomy. The Casimir energy is given 
by (13.151) for both the spin structures and it is precisely the same as (14.81) . The geometry 
of AdSs is believed to be stable under a' correction and the stability seems to continue 
to the standard orbifold. 

The thermal AdS orbifold has a fixed point at r = and there are closed strings 
in twisted sectors localized at the fixed point. With the periodic boundary condition 
for the fermions along the x-cycle, there are massless states localized at the fixed point, 
and the different vacua are dual to different excitations of these massless states. The 
excitation of localized massless states does not change the global geometry, thus all the 
geometries contribute equally to the partition function. In particular, the mass of the 
every geometries should be the same as ( 14.81) . which is also the same as the Casimir en- 
ergy for the every vacua. Therefore, we can say that the phase structure does not change 
even at the large 't Hooft coupling!^] For large k, the T-dual picture along the x-cycle 
is relevant, where the k NS5-branes are arranged in the dual x-cycle [33] ■ In the T-dual 



picture, the different excitations of massless states correspond to different configurations 
of k NS5-b r anes§ In particular the vacuum with the trivial holonomy corresponds to 
the configuration of k coincident NS5-branes. With the anti-periodic boundary condition, 
there are localized tachyons at the fixed point as in [U] . The different vacua are dual to dif- 
ferent condensations of these localized tachyons, which deform the geometry significantly 
from the thermal AdS orbifold. Since the configurations with tachyon are unstable, the 
relevant geometry should be the one without tachyon, which will be discussed in the next 
subsection. 

At high temperature, the dual geometry is the orbifold of the AdS-Schwarzschild black 
hole ([I~3D with ds| 3 replaced by ll£7J) . The mass of the black hole may be given by the 
expansion of the temperature as 

M = ~ ^ + 0{l) ' (4 ' 9) 



21 We are not sure whether the phase structure is the same in a middle value of the 't Hooft coupling, 
but we guess that this is indeed the case. 

22 It was shown in [S] that the dual geometries are also labelled by the integers (no, ■ • • , nk-i) at zero 
temperature. They discussed how to construct these geometries, where the NS5-branes are replaced by 
flux and the back- reaction of the flux is taken into account. 
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which is 3/4 times the energy given in (I3.3U|) as discussed in [2H]- Since the geometry does 
not have any fixed point, there are no light localized modes generically. Due to this fact, 
the orbifold of the black hole seems to be the most relevant one, even though the every 
vacua are degenerated at the zero 't Hooft coupling. At relatively small temperature or for 
large k, there may be nearly massless modes or tachyonic modes near the horizon. With 
the periodic boundary condition, the T-dual picture is more relevant for large k, where 
the shift of NS5-branes is given by the nearly massless modes. With the anti-periodic 
boundary condition, the orbifold of the AdS-Schwarzschild black hole may decay into a 
resolved AdS orbifold through the tachyon condensation. This type of geometry decay 
has been discussed in [3J1 [33] • According to them, a kind of black hole decays into a 
bubble of noting by a winding tachyon condensation. 



4.3 Localized tachyon condensation 

If the anti-periodic boundary condition is assigned for the fermions along the x-cycle 
then there are tachyonic modes at the fix point of the thermal AdS orbifold at low tem- 
perature. The geometry deformed by a tachyon condensation was proposed by [10J [11 



a0 

dr 2 T 2 

ds 2 = g{r)dt 2 + ^TjT^j + j[f(r){dx + cos^0) 2 + d9 2 + sin 2 6d<p 2 } , (4.10) 

with 

g(r) = r 2 + 1 , f(r) = 1 - £ , a 2 = - l) . (4.11) 

They call this geometry as Eguchi-Hanson soliton. We set k even in order to assign the 
anti-periodic boundary condition for the fermions and also k > 2 for a 2 > 0. Because of 
the boundary condition for the fermions, the non-trivial x-cycle can be pinched off. The 
boundary metric at r — ► oo is given by (14.71) with \ ~ X + 4:ir/k as supposed to be. The 
solution is regular everywhere including r = a, and there are no tachyonic modes. There 
is no geometry for r < a, and this region may be replaced by the tachyon state or the 
nothing state in the sense of [33] . 

The conserved mass of the deformed geometry was computed in [TUJ [TT] as 

M= 7r(3 - 4fl2) = - 7r(fc4 - 8fc2 + 4) (412) 
32A;G 5 128^5 ' K ' J 

where G5 = ir/(2N 2 ) in the gauge theory terminology. Since the non-trivial x-cycle is 
pinched off at r = a, the Wilson loop along the x-cycle can take non-trivial value according 
to the previous discussion on the Polyakov loop. This means that the symmetry is 
broken in this background, and this is consistent with the fact that the dual vacuum has 



23 See also 
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the trivial holonomy. It is amusing to notice that the mass of the Eguchi-Hanson soliton 
(I4.12p is about 3/4 times the Casimir energy of the gauge theory on the trivial holonomy 
vacuum (13.171) for large k. This reminds us of the famous 3/4 difference in the context of 
[30] mentioned above, but their origins are not directly related to each other. It is worth 
to study this issue furthermore. 

At enough high temperature, only the relevant geometry seems to be the orbifold of 
the AdS-Schwarzschild black hole. Therefore we may observe the phase transition between 
the Eguchi-Hanson soliton and the orbifold of the AdS-Schwarzschild black hole. As in 
the previous case the partition function may be obtained from the classical actions for 
these geometries in the classical approximation. The difference of these actions is given 
by 

V BH (r m ) - V EH (r m ) vr 2 r+(r j -r%- 2a 2 ) 
A 27rG 5 4fcG 5 (2r 2 +1) ' 1 } 

where the region of r < a is removed in the Eguchi-Hanson soliton. The critical temper- 
ature iS 

T c = 2 _±y^±K (4 14) 

V / Mu7rw) 

We can see from (I4.13P that the Eguchi-Hanson soliton is dominant at lower temperature 
and the orbifold of AdS-Schwarzschild black hole is dominant at higher temperature. 



5 Conclusion and discussions 

We have studied the thermodynamics of M = 4 U (N) super Yang-Mills theories on 
R x S^/Zfc with large N. The base manifold S^/Z^ has a non-trivial cycle along the 
X-direction of (14. 7p . and a non-trivial holonomy can be assigned along the non-trivial 
cycle. The theory has multi-vacua associated with the choice of holonomy labelled by k 
integer numbers (no, • ■ • , Uk-i)- We can assign the periodic and anti-periodic boundary 
conditions for the fermions along the non-trivial cycle. On a compact manifold, the Gauss 
constraint only allows gauge invariant operators, and due to this fact a phase transition 
occurs even in the zero 't Hooft coupling limit for large N. We have computed the partition 
function for the large N gauge theories with different holonomies by following the analysis 
in j2j [3], and examined the phase structure with a special care on the difference between 
the vacua. 

At low temperature, the most relevant contribution to the free energy comes form the 
Casimir energy. For the case with the periodic boundary condition, the Casimir energy 

24 Notice that the critical temperature of the Hawking-Page transition T c = 3/(2ir) is reproduced for 
a = 0. 
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does not depend on the choice of holonomy, and hence the vacua are degenerated. For 
the case with the anti-periodic boundary condition, the Casimir energy depends on the 
choice of holonomy, and the dominant contribution comes from the vacuum with the 
trivial holonomy V — 1. Near the critical temperature, we can obtain an approximate 
analytic expression of free energy by utilizing the Gross- Witten ansatz |17J. See Figure 
UJ For enough large k the case with the trivial holonomy seems to dominate for both the 
spin structures. At high temperature, the free energy is universal as in (13.301) . thus the 
vacua are degenerated for both the spin structures^ 

In the limit of large 't Hooft coupling, the dual gravity description is more appropriate. 
In the low temperature phase, the dual geometry is the orbifold of the thermal AdS space 
or its deformation. The orbifold of the thermal AdS space has the fixed point at r = 0, 
and there are closed string states localized at the fixed point. With the periodic boundary 
condition, the localized states are massless, and the excitation of these massless states 
leads to degenerated different geometries. In a T-dual picture, the different excitations 
correspond to different arrangements of k NS5-branes along the dual x-cycle. With the 
anti-periodic boundary condition, the localized states are tachyonic and a condensation of 
the tachyonic modes leads to the decay of geometry into the regularized geometry (14.101) 
called as the Eguchi-Hanson soliton [TUl ITT] . In the viewpoint of the dual gauge theory, 
the localized tachyon condensation is realized as the transition between different vacua. 
In particular, we have found that the mass of the Eguchi-Hanson soliton is about 3/4 
times the Casimir energy for the dual vacuum with the trivial holonomy. I would like to 
study the relation between the vacuum transition of the gauge theory and the RG-flow 
or time-dependent process among the geometries as in [§]. At enough high temperature, 
the dual geometry is the orbifold of the AdS-Schwarzschild black hole. Since there is no 
fixed point in this geometry, the orbifold seems to be the most relevant geometry. This 
implies that the phase structure would vary as the 't Hooft coupling is changed at high 
temperature. 

There are several theories similar to our orbifold gauge theories in the sense that the 
gauge theory has many vacua and its dual gravity description. One of them is (1+1) di- 
mensional large N gauge theories on a torus [371 138] , where non-trivial holonomy matrices 
can be assigned along the two cycles. For the case with the periodic boundary condition 
of the fermions along the spatial cycle, the eigenvalues of spatial holonomy matrix corre- 
spond to the positions of N DO-branes along the T-dual spatial cycle [371 EE]. I n the high 
temperature phase the eigenvalues are distributed uniformly, but in the low temperature 
phase the eigenvalues get together. This is related to the Gregory-Laflamme transition 



25 The masses of lightest twisted string modes are proportional to the horizon radius r+, and these 
modes become massless (or tachyonic) in the zero 't Hooft coupling limit since the radius behaves like 
r + ~ a' VAT '. From this reason we can assign non-trivial expectation values to the lightest modes and 
the vacua can be degenerated in the gauge theory description. 
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from a black string wrapped on a spatial cycle into a localized black hole . For the case 
with the anti-periodic boundary condition, the spatial cycle can shrink, and indeed the 
thermal AdSs is dominant in the low temperature phase. In this phase, the symmetry 
of the holonomy matrix along the spatial cycle has to be broken. In the high temperature 
phase, the BTZ black hole dominates, and the Z^v symmetry is preserved. The relation 
to our cases may be examined for large k limit, where the back-reaction of NS5-branes 
should be taken into account. I would like to study this relation, for instance, by taking 
a large N, k limit with keeping a ratio N/ k finite] 26 ! 

Other interesting theories are the plane wave matrix model [UJ and (2+1) dimensional 
super Yang-Mills theory on R x S 2 , which are obtained by truncating the J\f = 4 Yang- 
Mills theory on 1 x S 3 like our orbifold theory. These models share the same symmetry 
£77 (2 1 4) at zero temperature, and the gravity dual of these models was studied in [6]. The 
thermodynamics of plane wave matrix model has been studied in [4"2| |4"3"1 HH H5j 06], and, 
in particular, the different vacua was compared in [30] • The thermodynamics of large N 
gauge theory on R x S 2 should be also interesting. The relation among these models at 
zero temperature has been discussed in P, HZJ HEJ E] , and it is worth while investigating 
the relation among their thermodynamics. 

In this paper, we have taken the zero 't Hooft coupling limit A = 0, namely, the free 
theory limit, thus a next task is to study the effects of non-zero 't Hooft coupling. In the 
free theory limit, the phase transition is of the first order, but the inclusion of small A 
may change the order of phase transition as discussed in P, 09] . Moreover, we may be 
able to examine how the phase structure in the zero 't Hooft coupling limit continues to 
the one in the large 't Hooft coupling limit@ For example, the Casimir energy on the 
trivial vacuum with the anti-periodic boundary condition ( 13. 17ft is about 4/3 times the 
mass of the Eguchi-Hanson soliton (14.12p . thus one may wonder what would happen if we 
include A correction. In the high temperature phase, we have observed that the relevant 
geometry is the orbifold of the AdS-Schwarzschild black hole, which is dual to the Z fc 
symmetric vacuum. Since the vacua are degenerated in the zero 't Hooft coupling limit, 
the vacuum structure should depend on the 't Hooft coupling. We would like to study 
this issue as a future work. 
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A The partition function of single scalar particle 

As an example we compute here the partition function of single scalar particle. Let 
us define the following function on S 3 as 

f(x,y)=Trx E y 2 » . (A.l) 

Since a scalar is expanded by the scalar spherical harmonics Sj >m; rn(fl,), the above function 
can be easily evaluated as 

f(x, y) = x + 2x 2 (y + l/y) + 3x 3 (y 2 + 1 + 1/y 2 ) + ■■■ . (A.2) 

The orbifold case is given by restricting the modes into the ones invariant under the 
orbifold action (jOj) . 

In the case of the trivial holonomy V — 1, we have to project the modes into the ones 
with 2m = kZ. The projection leads to 



f(x 7 y ) = x + 3 x 3 + ... + Y / (y nk + l/y nk ) [ink + l)x nk+1 + (nk + 3)x nk+3 + ■••], (A.3) 

n=l 

where only terms proportional to y kZ are kept. Using a formula 

-(L-2) 
[1 - x 2 ) 2 



Lx L + (L + 2)x L+2 + ■■■= } L 2 f 2 )xL+2 , (A.4) 



the function can be written as 

I 3 oo / u , i \ nk+1 /„L _ 1 \™nfc+3 

/(*. v) = 0^ + Uv nt + Vt- r + ' (1 _ J )2 ' • (A.5) 

With Y^=i % n = - x) and Y%?=\ = X I0- ~ x ) 2 , we finally find for y = 1 

(x + x 3 )(l — x 2k ) + 2kx k (x — x 3 ) 

/(X ' 1)= (1- 1)2(1 • (A ' 6) 

This gives the expression in (12.51) with I = J. 

The case with non-trivial holonomy can be analyzed in a similar way. For a bi- 
fundamental scalar (rii, nj), the orbifold action is given by ( \2A\\ . In this case we keep the 
terms proportional to y kZ+L ^ and we find 



f(x, y) = J2 y nk+L [(nk + L + l)x nk+L+1 + (nk + L + 3)x nk+L+3 + • • ■ ] 

n=0 

(A.7) 

oo 

+ E y nk ~ L K nk ~ L + l)x nh ~ L+1 + (nk-L + 3)x nk ~ L+3 + ■■■]. 

n=l 

Making use of the formula ( 1A.4I) . we can reduce the above sum into the more simplified 
form in (12. 5p . 
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B An index for the supersymmetric orbifold theory 



In this appendix the index proposed in [19J is computed for Af = 4 U(N) super Yang- 
Mills theoryon K x S^/Z^ with the periodic boundary condition for the fermions along 
the x-cycleo The index in our case is defined as 

J = Ti(-l) F e-^t 2{E+n) v R ' 2 w R3 , S = E - 2ji - ^-R x - R 2 - ^R S • (B.l) 

Since only states satisfying H = contribute the index [19], we can set (3 — > oo. As before 
E is the energy and (Ji, j 2 ) are spins with respect to SU(2) 1 x SU (2) 2 isometry on S 3 . The 
symmetry of the theory is <S77(2|4) and R\, R 2 , R3 are .R-charges. With generic holonomy 
(no, • ■ • , Tik-x), the U(N) gauge symmetry is reduced to Y\j U(ni), and the states are in 
adjoint or bi-fundamental representation of the gauge group. As in [19, 20J or in (12.131) 
the index can be written by an integral of unitary matrices as 



2 = J[Y[dUj] exp 



00 1 

E E -f' J (t n , v\ ™ n ) T*(U?) MUT) 

I, J n=l 71 



(B.2) 



1 

where the coefficients f I,J are the indices for the single particles with representation 
(ni, fij). Notice that the index does not receive any corrections of 't Hooft coupling, since 
we are counting the modes protected by supersymmetry [19J. 

Let us compute the indices for single particles. We start from the scalar contribution. 
The theory includes three complex scalars X,Y,Z with (Ri, R 2 , R 2 ) = (0, 1, 0), (1, —1, 1), 
(1, 0, —1). The scalars are expanded by Sj im ,m(^), and the condition H = is satisfied by 
the modes Sjj/ 2jfh . The orbifold projection means 2m e L + fcZ for the L = J — I + A;Z 
(0 < L < k) sector, and the contribution can be computed as 

t 2 (t 3L + i^ k -^){v + l/w + w/v) 

(1 - t 6 )(l - t 3fc ) ' ( ' 

Since the gauge field does not have any .R-charges, the condition 5 = is satisfied by the 
modes V^j +1 y 2 lfh with 2m G L + kZ. Thus the contribution from the gauge field is 

t 6 (t 3L +t 3( fc -L) ) 

(1 - t 6 )(l - t 3fc ) ' 1 ' ' 

For the fermions, the condition H = can be satisfied by the modes F^, 2 m with 2m e 
L + kZ and (i?i, R 2 , R 2 ) = (1, —1, 0), (0, 1, —1), (0, 0, —1). The contribution from these 
fermions is 

t 4 (t 3L + t 3( - k -^)(l/v + w + v/w) 

(B.5) 



;i-t 6 )(l-t 3fc ) 



2S 



See also [52 [MIES]. 
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The other contribution comes from the modes F, 
(1,0,0) as 

5l.o 



j,(j-i)/2,m witn 2m G L+kZ and (R 1} R 2 , R 2 ) 

( f 3L +f 3(fc-L)) 



(B.6) 



(B.7) 



(1 - t 6 )(l - t 3& ) ' 
Interestingly the sum of all contributions can be factorized as 

£lJ , N r (t 3L + t 3 ( fc ^))(l-t 2 t;)(l-t 2 /w;)(l-t 2 w;/^ 

/■ (M,^) = ^0 (l-t«)(l-t3*) • 

For the indices to converge, we have to assign t 2 < 1, t 2 v < 1, t 2 /w < 1, t 2 w/v < 1, which 
implies 5jj — / 7 ' J (t, v , w) > 0. 

In order to perform the integral (IB.2I) . we assume that nj is very large or zero. Then we 
can replace the discrete eigenvalues by continuous ones with the densities p T (9i) satisfying 
J d9ip I (9i) — 1. In this term, the effective action is given by 



i,j 



n i 



V, 



n=l 



' ■ J = -(6i,j-f I ' J (t n ,v n ,w n )) 
n 



(B.8) 



where we denote the Fourier transform of p I {9j) as p n . The saddle point of the action is 
p ! n = for n 7^ as V^' J > 0, and the index is given by the determinant 

1 



rI,J\ 
'n i 



(B.9) 



„ det(Ei,jninjV n 

The determinant is complicated in general, but it could be written in a simpler form for 
several cases. 

One case is with the trivial holonomy V = 1. In this case, we have just lxl matrices, 
so the determinant is simply 

J = n (l-t 6 ")(l-t 3fc ") 

L } (l + t 3kn )(l-t 2n v n )(l-t 2n /w n ){l-t 2n w n /v n ) ' 1 ' ; 

The N dependent factor is removed by changing the normalization. Another interesting 
case may be with the symmetric holonomy nj = N/k for all I. In this case, it is useful 
to utilize a formula for a circulant determinant as in [20] 



fi f2 fz 

fk fl f2 
fk-l fk fl 



fk 
fk-l 
fk-2 

fl 



k-l 



U(h + ^ I f2 + ^ I h + --- + ^ k - 1)I fk ) 

1=0 



fB.lll 



f2 fz fi 

with uj = exp(27ri / k) . Using the identity 

/=0 L=0 



'1-t 



3k\ki 



1 - t 



6\k 



[i-t 3k y 



;b.i2) 
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the index is computed as 

1 ~ (1 _ t 2n v ny(l _ ^^n)^ _ iV/« B )' ' (B.13) 

We would like to compare these indices with the one from the gravity computation. 
However there is a subtle problem whether we should remove the contribution from the 
diagonal U(l) k ^ 1 part. In the case of AdS$ x S 5 /^, the U(l) k ~ 1 part of the dual Yang- 
Mills theory is removed to compare the indices with those from the gravity computation 
[20] . Suppose that the similar U(l) decoupling should be taken into account even in our 
case. The contribution from U(l) part is just the same as that of V — 1 case (IB.10I) . so 
after the subtraction we have 

/i j.3kn\2(-\ _i_ j.3kn\k— 1 

j _ tt \± 1 ) \ L X i CB 141 

A n A (1 -t 6n ) fe - 1 (l -<3*n)*-l(l -t 2n v n )(l -t 2n /w n ){l -t 2n W n /v n ) ' 1 ' ' 

We would like to study the U(l) problem furthermore to examine whether this is indeed 
the case. 

The dual geometry is the orbifold AdS^/Zk x S 5 with the excitation of localized string 
states at the fixed point, and the different excitations correspond to different vacua of the 
gauge theory. The index for the supergravity on AdS^fL^ x S* 5 can be computed by acting 
the orbifold projection to the index in the covering space case P3EEII]- The single-particle 
index can be computed as 

k r t 2 v t 2 /w t 2 w/v t 3 ^ 1 tf/u 1 



K 1=1 L 



+ z zzn— + 



t 2 v l-t 2 /w l-t 2 w/v 1-tV l-tf/uj 1 
t 2 v t 2 /w t 2 w/v 2t 3k 



;b.is) 



l-t 2 v l-t 2 /w l-t 2 w/v l-t 3fc ' 
and the index including the multi-particle contribution is given as 

1 = 6XP S n 1 '^'^' ^ = S (1 - P^Afl - tW w n Kl _ t 2n w n /vn) " ( B ' 16 ) 



n=l n=l 



The index is a slightly different form the both of (IB.lOj) and (1B.14j) . and the difference 



should be interpreted as the contribution from the localized string states. It is desirable 
to include these stringy contributions to the index from the AdS side and compare it with 
the index of the gauge theory on the corresponding vacuum. 
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